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We present a general formalism for describing stimulated Raman adiabatic passage in a multi- 
level atom. The atom is assumed to have two ground state manifolds a and b and an excited state 
manifold e, and the adiabatic passage is carried out by resonantly driving the a — e and b — e 
transitions with time-dependent fields. Our formalism gives a complete description of the adiabatic 
passage process, and can be applied to systems with arbitrary numbers of degenerate states in each 
manifold and arbitrary couplings of the a — e and b — e transitions. We illustrate the formalism by 
00 ' applying it to both a simple toy model and to adiabatic passage in the Cesium atom. 
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A number of techniques for coherently manipulating atomic systems are based on the idea of adiabatically varying 
the parameters of a time-dependent Hamiltonian. One such example is the technique of stimulated Raman adiabatic 
passage (STIRAP), which can be used to coherently transfer population between the two ground states of a three- level 
i-M ' atom in the lambda configuration [EH 0]. The STIRAP technique involves dr iving the atom with a pair of time- 
dependent fields that couple the two ground states to the excited state, and relies on the fact that the Hamiltonian 
for the system has an instantaneous eigenstate, called the dark state, that contains no excited state component. In 
general this dark state is a superposition of the two ground states, where the overlap of the dark state with each 
ground state depends on the ratio of the powers in the two fields. Population is transfered from the first ground state 
to the second by adiabatically varying these powers in such a way that initially the dark state overlaps entirely with 
'— the first ground state, during the process the dark state is a superposition of the two ground states, and after the 
process has been completed the dark state overla ps e ntirely with the second ground state. The STIRAP technique 
has been studied experimentally [3) S S 0, B B lid ] , and has formed the basis for numerous theoretical proposals 
■ [HI [H, [H, [lij. In addition, STIRAP has been generalized to multi- level systems [ljl [l?], EH], opening up new 
Q\ | possibilities for creating superposition states [H, 0, HH, [H, HH . 

t — ■ Here we present a general formalism for describing STIRAP in a multi-level atom. The atom is assumed to have two 
ground state manifolds a and b and an excited state manifold e, where each manifold consists of an arbitrary number 
. of degenerate states, and the adiabatic passage is carried out by resonantly driving the a — e and b — e transitions 
C^) " with time-dependent fields. The Hamiltonian for the system is 

oo ! 

H=-cos0(A b +Al)- - sind (A a +Al), (1) 



where Q is a time-independent Rabi frequency that characterizes the strength of the fields, 9 is a time-dependent 
parameter that is varied to carry out the adiabatic passage, and A a and Ab are lowering operators that connect states 
in manifold e to states in manifolds a and b. We sweep 9 from to tt/2 to transfer states from manifold a to manifold 
b, and we sweep 9 from 7r/2 to to transfer states from manifold b to manifold a. 

Let us denote the Hilbert spaces for manifolds a, b, and e by 7i a , Tib-, and TL e . We will show that TL a and TLb can 
be decomposed as follows: 

n a = n x a ®n d a ®K (2) 
n b = nt®nt®n' b . (3) 

Under adiabatic passage from a to b, states in Tt* coherently evolve into states in states in 7if remain unchanged, 
and states in TL' a are driven to the excited state manifold and scatter incoherently. Similarly, for adiabatic passage 
from b to a, states in H$ coherently evolve into states in H^, states in remain unchanged, and states in Ti' b are 
driven to the excited state manifold and scatter incoherently. The coherent evolution from to H.^ is described 
by a unitary transformation U : — * H^, and the coherent evolution from H£ to is described by the unitary 
transformation U' : TL^ — > Ti-a- Given arbitrary operators A a and Ab, our goal is to perform the Hilbert space 
decompositions described in equations @ and ([3]) and to calculate the unitary transformation U. 

The paper is organized as follows. In section[TTl we show how to preform the Hilbert space decompositions described 
in equations @ and ([3]). In section HTT1 we calculate the unitary transformation U : Hfc, and consider in detail 
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the special case dimH^ = dim7Y^ = 2. The formalism developed in sections HT1 and HTT1 gives a complete description 
of the adiabatic transfer process, and in section ITVl we illustrate this formalism by applying it to a simple toy model. 
Finally, in section |V] we use the formalism to analyze adiabatic passage in the Cesium atom. 

II. HILBERT SPACE DECOMPOSITIONS 

We will first show how to perform the Hilbert space decompositions described in equations ^ and © for a given 
pair of atomic lowering operators A a and Aj,. Let H d be the space of states in manifold a that are dark to the a — ► e 
transition, let Hf be the space of states in manifold b that are dark to the b — > e transition, and let Ttf be the space 
of states in manifold e that are dark to the e — > b transition. Note that H d , H d , and Hf are just the null spaces of 
the operators A a , A\, and Ab: 

K = {\tpa) S H a | Al\lp a ) = 0} (4) 

Hf = {\i) b } e H b I = 0} (5) 

ttf = {IV'e) e H e I A b |^ e ) = 0}. (6) 

Define H^, H^, and Hjr to be the complements of these spaces in H a , Hb, and H e : 

H a = H^® H d a (7) 
H b = Hi®H d b (8) 
W e = Hi®H d e . (9) 

Clearly, states in and Hf are dark states of H . In addition, H has dark states of the form 

| A) = cos0|Va> + sin <f>\i/j b }, (10) 

where \ip a ) & H^ and \ipb) € We will call such states lambda dark states, and say that \tp a ) and \tpb) form a 

lambda pair. Since i?|A) = 0, we have that 

cos 9 sin <$>A\ \ ipb) = sin#cos0 A a \ip a )- (11) 

This equation must hold for all values of 9, so <f> must be related to 9 by 

tan cj) = A tan 9 (12) 

for some value A. Thus, states \ip a ) and \ipb) are related by 

\Al\ih) = A\\i> a ). (13) 

Note that by using equation (fl"2")) . we can also express |A) as 

| A) = (cos 2 9 + A 2 sin 2 <9)- 1/2 (cos6> \tp a ) + Xsin9 \tp h ))- (14) 

Define H^ to be the space of states \tp a ) € such that A a \ip a ) € 7i^-, and define 7^ to be the complement of H* 
inHt- 

Hi = H x a ®H' a . (15) 

Note that from the definitions of H' a and Hf., it follows that dim 7^ < dimH^. 

We will now show that for every normalized state \ip a ) G H x , we can construct a normalized state \ijjb) G 
such that |f/> a ) and j^b) form a lambda pair. First, note that we can view the raising operator A\ as a mapping 
A\ : H a — * He ■ Since the image of H x under A\ lies entirely in H^r , we can define a new mapping i?| : ft* — > by 

Rl\4> a ) = 4|^a), (16) 

where |i/> a ) is an arbitrary state in 7i^. Similarly, given the definitions of H^ and Ti^i we can define mappings 
Rb : Hi — ► We 1 and R b : H^ ^ H^ by 

i?t|V- b ) = At|V%) (17) 

ifelVe) = Afcl^e), (18) 
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where \4>b) and \ip e ) are arbitrary states in 7i b and H^. Note that the only state in that lies in the null space of 
A\\ is the null state, the only state in H. b that lies in the null space of A\ is the null state, and the only state in Tijr 
that lies in null space of At, is the null state; thus, the null spaces of R\, R b , and Rb are trivial. This means that the 
mapping R) b Rb : TLj: — > TLj: is invertible, and we can define a mapping M : — > by 

M = RbiRlRb)- 1 Rl (19) 

Define to be the image of M, and define 7^ to be the complement of Ti b in Ti b : 

Hi =H£® H' b . (20) 

Note that because the null spaces of Rl, Rl, and R b are all trivial, the null space of M is also trivial, and therefore 
dimH b = dimH^. Given a normalized state \ip a ) S H*, define A by 

A= {iP^M^a) 1 ' 2 , (21) 

and define a normalized state \ipb) G Ti b by 

Hb) = \M\1> a ). (22) 

Using equations (fH)|) and (f2"2"|) . we find that 

A4|^) = AflJ|^) = R\M\^ a ) = R\R b (RlRb)- 1 R\\i>a) = R^a) = 4lV>a>- (23) 

Thus, states \4>b) and \ip a ) satisfy equation (TT3")) and therefore form a lambda pair. 

We claim that none of the states in T-i' a can form lambda pairs with states in H. b . To see this, consider a superposition 
of states \tp a ) £ H' a and \ipb) € H b : 

\4>)=Ca\^a)+Cb\^b)- (24) 

Because of the way we have defined H' a , there must be a state \ip e ) € Tif such that 

(V> e |4 \^ a ) ? 0. (25) 



Since |?/>e) G Tig we have that Aftj^e) = 0, so 



(Ve|4lV>b) - 0, (26) 



and therefore 



(i> e \H\<p) = ~ sin9 W e \Al\i> a )c a . (27) 



Thus, for c a 7^ and sin 6* ^ we have that H\4>) ^ 0, so \<f>) cannot be a lambda dark state. 

Since all the states in form lambda pairs with states in Tt b , and none of the states in H' a form lambda pairs 
with states in H b , it follows that none of the states in H! b form lambda pairs with states in H^. Thus, states in H' a 
and TL b do not form dark states of H , and under adiabatic passage they are driven to the excited state manifold and 
scatter incoherently. 



III. UNITARY TRANSFORMATION 



In the previous section we defined Hilbert spaces and Ti. b , which are subspaces of the total Hilbert spaces for 
ground state manifolds a and b. We will now show that under adiabatic passage from a to 6 states in H.^ coherently 
evolve into states in H b , and we will derive the unitary transformation U : — > Ti. b that describes this evolution. 

First, choose an orthonormal basis {l^ai), • • • , for 7i^, and use the mapping M given in equation (fT9|) to 

construct a normalized basis {\ipbi}, ■ • • > \ipbn)} f° r H b , where \4>bj) is defined by 

lV>w) = ^-Mtyaj) (28) 
A, = (ip a j\M^ Mlipaj) 1 / 2 . (29) 



Note that the basis states for Tt£ will not necessarily be mutually orthogonal. The states \t^ a j) and \i/>bj) form a 
lambda pair, and define a lambda dark state that is given by 

\A,{9)} = cos^(0) \il> aj ) + smfaid) |^-), (30) 

where 

(j>j(6) = tan" 1 (A.; tan 0). (31) 

The states { | Ai) , • • • , |A„)} form a basis for the lambda dark states, so we can express a general lambda dark state 
| A) as 



\A) = ^2c j \A j ). (32) 



In the adiabatic limit, the time evolution of |A) is given by Schrodinger's equation: 

i^|A) = H|A) = 0. (33) 

Thus, substituting for |A), we find that 

J2(c j \A j )+c j ^ j \A j )) = Q, (34) 

3 

where the dots indicate derivatives with respect to 9, and where we have defined 

Taking the inner product of equation (|34|) with ( A& | , we obtain 

^(cj (A fe |A,-) + cj fa (A fc |A,)) = 0. (36) 

3 

Since {|Ai), • • • , |A„)} is a complete basis for the lambda dark states, the matrix (A/-|Aj) is invertible; we will denote 
its inverse by L^: 

^L ik {\ k \\ j )=8 ij . (37) 
k 

If we multiply equation (|36p by La- and then sum over k, we obtain the following equations of motion for the expansion 
coefficients a: 

di = -^2^2Uk (A k \Aj) <pj cj. (38) 

3 k 

Note that we can express an arbitrary state \ip a ) G W„ as 

\i>a) = ^>(0) IVW> = E C J(°) l A ^°))- ( 39 ) 
3 3 

for some set of amplitudes Cj(Q). We can then integrate equation (|38|) subject to these initial conditions to obtain a 
state 

U\4>a) = 5>(tt/2) |A j -(tt/2)> = X>(tt/2) (40) 

3 3 

This defines the unitary transformation U. 
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A. Example: dim'W* = dim 7^ = 2 

As an example, we will write down the equations of motion explicitly for the case dim7i^ = dimWt = 2. We choose 
an orthonormal basis of states {IVvi), IVV2)} for use equation (|2"5jl to obtain a basis of states {\ipbi}, IVte)} for 
H.£, and use equation (|2"5|) to obtain the values Ai and A2. The matrices (A^Aj) and (Aj,|Aj) are given by 

(A k \Ai) = * sin< ^ sin M (41) 

N 1 31 \ z* sin 4>i sm 4> 2 I J v ; 

(A fc |A,) = ( . ,° . , * sin ^ cos M (42 ) 
\ 2* cos 0i sin 2 / v ; 



where 



The inverse of (Afc|Aj) is 



= (->Pu\ipb2). (43) 



r /-, 1 12 • 2 / • 2 , \-i / 1 —z sin 0i sin (/> 2 \ / a a\ 

= (l-k| 2 sin 2 lS m 2 2 ) \_ z * shl(j)ishl(j)2 ? 92 J- (44) 

Substituting these matrices into equation f|38|) , we find that the equations of motion for c\ and C2 are 

6 }\ =(l-|z| 2 sin 2 1 sin 2 2 )- 1 p z|2sin / lC °^l si l l2 ^ . l2 " Z f 01 T S02 2, 1 ( V (45) 
c 2 / v 1 1 y V cos 0i sm 02 |z| i sin0 2 cos02sm 2 0i J \ 2 c 2 / 

and from equation (|31[) we have that 

cos fe = (cos 2 9 + Xl sin 2 6»)~ 1/2 cos (9 (46) 

sin0 fe = (cos 2 (9 + A^sin 2 6»)- 1/2 A fc sin6l (47) 

fc = A fc (cos 2 6» + A^sin 2 6»)" 1 . (48) 

Given Ai, A2, and z, these equations of motion can be integrated to obtain the unitary transformation U : — > W^. 
For certain special cases, we can perform the integration analytically and write down the explicit form of U . 

First, suppose \ipbi) an d I^m) are orthogonal. Then z = and the equations of motion reduce to b\ — c 2 — 0, so 
the unitary transformation U : Ti^ — * H-^ is given by 

Utyal) = \M (49) 
U\i> a2 ) = \tp b 2). (50) 

Next, suppose Ai = A 2 , so there is a single angle = 0i = 02 that characterizes both lambda dark states. Also, 
for simplicity, assume that z is real. Then we can express the equations of motion as 

4)=^w(i-.w«->( 2 -;< 1 ^)(;;), (5i) 

where primes denote derivatives with respect to 0. We can decouple these equations by defining new variables 
r)± = ci ± c 2 : 

rj' ± = =Fzsin0cos0(l ± z sin 2 0) _1 rj±. (52) 
These equations can be integrated to give 

77±(0) = (l±z sin 2 0)- 1 / 2 7 ]± (O). (53) 

Thus, 

ci(0) = i[(l + zsin 2 0)- 1 / 2 ( Cl (O) + c 2 (O)) + (l-zsin 2 0)- 1 / 2 ( Cl (O)-c2(O))] (54) 
c 2 (0) = i[(l + zsin 2 0)- 1 / 2 ( Cl (O) + C2(O))-(l-zsin 2 0)- 1 / 2 ( Cl (O)-c2(O))]. (55) 
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FIG. 1: Level diagram for the example system. Horizontal lines indicate states in ground state manifolds a, b and excited state 
manifold e; diagonal lines indicate transitions coupled by A a and A\\ vertical lines indicate transitions coupled by Ab and A b . 

For adiabatic passage from a to 6 we sweep from to 7r/2, and from equation (|3Tj) we see that c/> also sweeps from 
to 7r/2. Thus, after the adiabatic passage has been completed 4> = ir/2, and the unitary transformation U : — > H^ 
is given by 

U\ip al ) = a\i> b i)+p\il! b2 ) (56) 
U\ip a2 ) = 0\^bi)+a\^ b2 ), (57) 

where 



-l/2 + (1 _ z) -l/2) 



(58) 
(59) 



IV. EXAMPLE SYSTEM 



We will now illustrate the formalism developed in sections HT1 and ITU with a simple example. Consider an atom that 
has the following internal states: 



H a = span{|a + ), |a ), |a_)} 
H b = span{|6+), |6o), |6_)} 
H e = span{|e+), |e ), |e_)}. 



(60) 
(61) 
(62) 



Define atomic lowering operators A a and A b by 



A b = |6 + )(e + | + |6 )<eo| + |6-)<e-| 

A a = |a + )«e + | + (e |) + |a_)((e_| + ( eo |). 



(63) 
(64) 



The transitions coupled by these operators are shown in Figure [T] 

We will first apply the results of section |TT] to find the decompositions of Hilbert spaces TC a and Hb- There is a 
single dark state |ao) for the a — > e transition, and there are no dark states for the b — > e and e — > b transitions, so 



K 

njd 

n b 
Hi 

H d e 

Hi 



span{|a )} 
span{|a±)} 

{} 
Hb 

{} 

H e - 



Since A\\tj} a ) € TLjr for every state \tp a ) & Hi, we find that 



H' 



Hi 
{}■ 



span{|a±)} 



(65) 
(66) 
(67) 
(68) 
(69) 
(70) 



(71) 
(72) 
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The mapping M : — > Ti.^ is given by 

M = (|6+) + |6o))(o+| + (|6_) + |6o))<o_|. (73) 
We can use this mapping to define states \B±) that form lambda pairs with states \a±): 

\B±) = ^M\a ± ) = (1/V2)(|fo ) + |&±», (74) 

where 

A± = (a±|M t M\a ± y/ 2 = y/2. (75) 
Note that \B+) and \BJ) are not orthogonal: 

z = {B + \B_) = 1/2. (76) 

From the results of section [Til it follows that 

|A±) = cos <P(6) \a±) + sin <P(d)\B±) (77) 

are lambda dark states of H, where 

(l>(9) =tan" 1 (\/2tan6l). (78) 

Define a state \s) in Tib that is orthogonal to both \B + ) and 

| S ) = (l/V3)(|6_>-|6 ) + |6 + ». (79) 

Since the image of M is span{|£?±)}, we find that 

H x h = span{|S ± )} (80) 
n' b = span{| S )}. (81) 



Now that we have decomposed the Hilbert spaces Tt a and Tib, let us apply the results of section HTT1 to calculate the 

A 

b ■ 



unitary transformation U : —> H^. We first note that a general lambda dark state of H is a superposition of |A 



and |A_): 

|A)=ci|A+)+c 2 |A_). (82) 

Since A + = A_, the equations of motion for C\ and C2 are given by equation (|5ip . As was shown in section fill Al these 
equations of motion can be integrated to give U : — > H £ : 

U\a+) = a\B + )+f3\B_) (83) 
U\a-) = P\B+)+a\B_), (84) 

where a and /3 are given by equations (|58p and with z = 1/2. 

The Hilbert space decompositions and unitary transformation U give a complete description of the behavior of the 
example system under adiabatic passage: for adiabatic passage from a to b, states \a±) coherently evolve into states 
U\a±) and state |ao) remains unchanged; for adiabatic passage from b to a, states U\a±) coherently evolve into states 
\a±), and state \s) is driven to the excited state manifold and scatters incoherently. 

It is interesting to note that if we start in state |a+) and perform adiabatic passage from a to b, population is 
transferred to |6_) despite the fact that (b_\Ab A\\a + ) — 0. The population reaches state |6_) by passing through 
state |ffl-); at time t the population in |a_) is given by 

P(|a_)) = |c 2 (0)| 2 cos 2 0, (85) 
where C2(<j>) is given by equation ([54|l with ci(0) = 1, C2(0) = 0, and where (f> — tan _1 (v / 2 tan 6(t)). 
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FIG. 2: Level diagram for Cesium. The two hyperfine ground state manifolds QSi/2,F = 3 and QSi/2,F — 4 correspond 
to manifolds a and b in our theory, and one of the excited state hyperfine manifolds, either Q>P 3 /2,F' = 3 or 6P 3 / 2 ,F' — 4, 
corresponds to manifold e. 



V. ADIABATIC PASSAGE IN CESIUM 



In order to show how the formalism developed in sections |TT] and IIIII applies in a more physical context, we will 
use it to analyze adiabatic passage in the Cesium atom. A level diagram for Cesium is shown in Figure [21 the two 
hyperfine ground state manifolds 6S , 1 / 2 ,-F = 3 and QSx/2,F — 4 correspond to manifolds a and 6, and one of the 
excited state hyperfine manifolds, cither 6P 3 / 2 ,F' — 3 or 6P 3 / 2 ,F' — 4, corresponds to manifold e. For simplicity, 
we will denote the 6Si/2,F = 3 and 6S±/2,F = 4 ground state manifolds by g% and g^, and the 6P^/2,F' = 3 and 
6P3/2, F' = 4 excited state manifolds by e 3 and 64. 

Adiabatic passage between manifolds (73 and 54 is performed by driving the atom with a pair of classical fields that 
connect these manifolds to the excited state manifold e.p>- The field driving the <-> ep> transition corresponds to 
an atomic lowering operator that is given by 



A F = A 



F ■ £F, 



(86) 



where If is the polarization of the field. The operator Ap is defined by 



4 = EE E P(F',F){F',m , \l,q;F,m)\F , ,m'){F,m\e* q 

m' m q— — l 



(87) 



where (F 1 , m'\l, q; F, m) is the Clebsch-Gordan coefficient that connects ground state \F,m) to excited state \F',m') 
via polarization e*, 

e±i = T^(x±iy) 
e = z, 

is a orthonormal basis of polarization vectors, and f3(F' , F) is given by the following table: 



(89) 



F' 


3 


3 4 4 


F 


3 


4 3 4 


(3(F\F) 


\/374 


s/T/l y/5/12 y/7/12 



(90) 



In the following sections we apply the results of section |TT] to perform the Hilbert space decompositions given in 
equations @ and for adiabatic passage via both the F' = 3 and F' — 4 excited state manifolds. 
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A. Adiabatic passage via the F' = 3 manifold 

For the g 4 <-» e 3 transition there are two dark states |c?4 e ) , |<i 4o ) in g 4 (these dark states are calculated in Appendix 
IA 2[) , and there are no dark states in e 3 , so 

Hi = span{|d 4e ), \d 4o )} (91) 
H d e = {}. (92) 

For the g 3 <-> e 3 transition there is a single dark state ^3) in 53 (this dark state is calculated in Appendix I A ip , so 

Hi = span{|d 3 )}. (93) 

From the dimensions of Hf , Hf, and , we can work out the dimensions of all the Hilbert spaces in the decompositions 
of H3 and H4. First, note that because dirnHf = and dimH 3 < dimHf, we must have that dim7Y 3 = 0. Next, 
since the dimensions of H 3 , H%, and H 3 must sum to 7, we find that dimHg = 6. Finally, since dimH 4 = dimWg = 6 
and the dimensions of H\ , Hf, and H 4 must sum to 9, we find that dimH 4 = 1. Thus, the dimensions of the Hilbert 
spaces are 



F 


dimW£ 


dimH d F 


dhnH'p 


3 


6 


1 





4 


6 


2 


1 



Note that these dimensions are independent of the polarizations of the classical fields. 



1. Example: £3 = £4 

The Hilbert space decompositions are 



Ht 


= span{|3,0)} 






= {} 




H3 


= span{3,±l), 


3,±2), |3,±3)} 


n/d 
/T 4 


= span{|4,±4)} 




7^4 


= span{|4,0)} 




njX 

rl 4 


= span{|4,±l), 


4, ±2), |4,±3)} 



(94) 
(95) 
(96) 
(97) 
(98) 
(99) 



States |3,m) and |4, m) are lambda pairs for m = ±1,±2,±3. Since the lambda dark states corresponding to these 
pairs are mutually orthogonal, equation (|3"8"]) decouples and the unitary transformation U : H 3 —*■ H\ can be written 
down explicitly: 

C/ = ^|4,m)(3,m|, (100) 

m 

where the sum is taken over m = ±1, ±2, ±3. 



2. Example: £3 = z, £4 = x 

The Hilbert space decompositions are 



Hi = span{|3,0)} (101) 

Hi = {} (102) 

H$ = span{|3,±l), |3,±2), |3,±3)} (103) 

Hf = span{|d 4e ), |d 4o )} (104) 

H' 4 = span{|s)} (105) 

H4 = span{|B ±1 ), \B ±2 ), \B ±3 )}, (106) 
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where 

| S ) = (1/n/2)(|4,1)-|4,-1)) (107) 

and the dark states are given by 

\d ie ) = (l/8\/2)(|4,4) + VM|4,2) + \/70|4,0)+'v / 28|4,-2) + |4,-4)) (108) 
\d 4o ) = (1/4) (|4, 3) + V7 14, 1) + V7 j4, -1) + [4, -3>). (109) 

We have defined states |B±i), \B±2), and \B±s) that form lambda pairs with states |3, ±1), |3,±2), and |3,±3); they 
are given by 

|B± 3 ) = (l/8V254)(127|4,±4)-2V7|4,±2)~y70|4,0)-2V7|4,T2)-|4,T4)) (110) 
\B ±2 ) = (l/4V^5)(15|4,±3)->/7|4,±l)->/7|4,Tl>-|4,T3» (111) 
\B±i) = (1/24V^)(99|4,±4) + 198v^|4,±2)-29v / 70|4 ! 0)-58v / 7|4,=f2)-29|4,=f4)), (112) 

and the corresponding A values are 



A ±3 = v /!0287/1792 (113) 

A ±2 = ^45/14 (114) 

A ±1 = v/8613/8960. (115) 

The unitary transformation U : H% — > Ti\ can be obtained by numerically integrating equation (|38p using the above 
expressions for the states and A values. We can say something about the structure of this unitary transformation by 
noting that the polarizations of the classical fields impose the selection rule Am = ±1. Thus, it is useful to decompose 
H% and H* into subspaces of even and odd Zeeman states: 

Hi = H$ e ®H%° (116) 

n\ = n\ e ®n\°, (117) 

where 

U\ e = span{|3,±2)} (118) 

U\° = span{|3,±l), |3,±3>} (119) 

U\° = span{|B ±2 >} (120) 

U\ e = span{|B ±1 ), |B ±3 )}. (121) 

The selection rule implies that U map s H^ e to H\° and H^° to U\ e . Since dimH^ = dim 7^° = 2 and A +2 = A_ 2 , 

we can use the results of section MI Al to write down the unitary transformation U : H% e — > H\ ° explicitly: 

U\3,+2) = a\B +2 ) + p\B- 2 ) (122) 

U\3,-2) = p\B +2 ) + a\B_ 2 ), (123) 

where a and (5 are given by equations (|55|) and (j59")l with z = {B +2 \B_ 2 ) = —1/15. 

B. Adiabatic passage via the F' = 4 manifold 

For the 54 *-* e 4 transition there is one dark state ^4) in 54 and one dark state \d e ) in (these dark states are 
calculated in Appendix I A ip . so 

H% = span{|d 4 )} (124) 

Hi = span{|4)}- (125) 

For the 53 <-> e 4 transition there are no dark states in (73, so 

H d 3 = {}. (126) 
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As with adiabatic passage via F' = 3, we can use the dimensions of H 3 , Ttf, and H d to say something about the 
dimensions of the other Hilbcrt spaces in the decomposition of H 3 and H^. For adiabatic passage via F' = 4, however, 
there are two separate cases to consider: since dim Tig < dim7if and dimH d — 1, we find that dimH 3 can be either 
or 1. For dim7i 3 = the dimensions of the Hilbert spaces are 



F 


dimW^ dimH d F dimH' F 


3 
4 


7 
7 1 1 



and for dim H' 3 — 1 the dimensions of the Hilbert spaces are 



F 


dim«£ 


dirnHj? 


dim7Y^ 


3 


6 





1 


4 


6 


1 


2 



Thus, for adiabatic passage via F' = 4 the dimensions of the Hilbert spaces depend on the polarizations of the classical 
fields. 



1. Example: £3 = £4 = I 
For this example dim7i 3 = 1, and the Hilbert space decompositions are 





= {} 


(127) 


K 


= span{3,0)} 


(128) 




= span{|3,±l), |3,±2), |3,±3)} 


(129) 


n_/d. 


= span{|4,0)} 


(130) 


K 


= span{4,±4)} 


(131) 


n_/X 
/T4 


= span{|4,±l), |4,±2), |4,±3)}. 


(132) 



States |3,m) and |4, m) are lambda pairs for m = ±1,±2,±3. Since the lambda dark states corresponding to these 
pairs are mutually orthogonal, equation (|38j) decouples, and the unitary transformation U : H. 3 — > H.* can be written 
down explicitly: 

(7 = ^|4,m)(3,m|, (133) 

m 

where the sum is taken over m — ±1,±2,±3. 



2. Example: £3 = x, £4 = z 



For this example dim7i 3 = 1, and the Hilbert space decompositions are 



where 



and 



H d 3 


= {} 






= span{|->} 




H 3 


= span{|3, ±2) 


|3,±3), |3,0), |+)} 


H d 4 


= span{ 4, 0)} 




H\ 


= span{|s e ), \s 




/T 4 


= span{|B ±2 ), 


|B ±3 ), \B ), \B + )}, 



|±) = (1/^)(|3,1) ±|3,-1)) 

\s e ) = (l/4)(|4,4) + V7|4,2)-V7|4,-2)-|4,-4)) 

|s ) = (l/V^)(3|4,3)+V7|4,l)->/7|4,-l)-3|4,-3)). 



(134) 
(135) 
(136) 
(137) 
(138) 
(139) 

(140) 

(141) 
(142) 
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We have defined states \B+) 
1 3, ±3); they are given by 



\Bq), \B±2), and \B± 3 ) that form lambda pairs with states 



\B±s) 
\B± 2 ) 
\Bo) 
\B + ) 



(1/V8) (V7|4,±4)-|4,±2)) 
(1/4)(V7|4,±3)-3|4,±1)) 
(1/V2)(|4,1) + |4,-1)) 
(l/V2)(|4,2) + |4,-2)), 



and the corresponding A values are 



A ±3 - v/25/49 



A ±2 = v/200/147 
A = V25CV49 
A+ = V375/392. 

Ti.^ can be obtained by numerically integrating equation 



|3,0), |3,±2), and 

(143) 
(144) 
(145) 
(146) 

(147) 
(148) 
(149) 
(150) 



The unitary transformation U : 7i 3 — ► TL^ can be obtained by numerically integrating equation (|38p using the above 
expressions for the states and A values. As with the example given in section IV A 21 the polarizations of the classical 
fields impose the selection rule Am = ±1, so it is useful to decompose TL 3 and 7i\ into subspaces of even and odd 
Zeeman states: 



At 4 



= n 



3 

Ae 



where 



n 
n 
n 



Xe 



\o 



\o 



The selection rule implies that U maps H3 



Al 4 = 

to ni° 



span{|3,0), |3,±2)} 
span{|+), |3,±3)} 
spanjlBo), \B ±2 )} 
span{|B+), \B ±3 )}. 

and H§° to n\ e . 



(151) 
(152) 

(153) 
(154) 
(155) 
(156) 



3. Example: £3 = £4 = e+i 



For this example dim?^ = 0, and the Hilbert space decompositions are 

nl = {} 

H' 3 = span{} 



/T4 
n_/\ 

Ha 



span{|3,0), |3,±1), |3,±2), |3,±3)} 

span{|4,4>} 

span{|4, -4)} 

span{|4,0), |4,±1), |4,±2), |4,±3)}. 



(157) 
(158) 
(159) 
(160) 
(161) 
(162) 



States |3, m) and |4, m) are lambda pairs for m — 0, ±1, ±2, ±3. Since the lambda dark states corresponding to these 
pairs are mutual 
down explicitly: 



pairs are mutually orthogonal, equation (|38|) decouples, and the unitary transformation U : H.3 — can be written 



U 



^2 |4,m)(3,m|. 



(163) 



m=— 3 



VI. CONCLUSION 



We have presented a general formalism for describing stimulated Raman adiabatic passage in a multi-level atom 
with ground state manifolds a and b and excited state manifold e, where each manifold consists of an arbitrary number 
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of degenerate states. This formalism describes how the Hilbert spaces for manifolds a and b decompose into subspaces, 
each of which evolves in its own characteristic way under adiabatic passage. In particular, we identified subspaces TL^ 
and Tib that coherently evolve into one another under adiabatic passage, and calculated the unitary transformation 
that describes this evolution. We have illustrated the formalism with a simple toy model, and used it to analyze 
adiabatic passage in the Cesium atom. The formalism gives a complete description of the adiabatic passage process, 
and should be useful for analyzing adiabatic passage in a wide variety of atomic systems. 
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APPENDIX A: DARK STATES 

Here we calculate the dark states for an optical field driving anF->Forf->F-l transition. We will choose a 
coordinate system such that the field propagates along the z axis, and define a basis of Zeeman states relative to this 
quantization axis. The polarization e of the field can be expressed as 

e = a+e+i + a_e_i, (Al) 

where e±\ are given by equation (I88|) . If \d) is a dark state of the field, then 

e-I t |d)=0. (A2) 

Thus, for every state \F' , m!) in the excited state manifold, we have that 

J2(F',m'\e-l^F,m)c m = 0, (A3) 

rn 

where 

c m = (F,m\d). (A4) 
Using equation (|87|) to substitute for A F , we obtain the following dark state conditions: 

{F',m\l,+l;F,m- l)a+c m _i + {F',m\l, -1; F, m + 1) ct-c m+1 = (A5) 
for —F' < m < F'. We will now consider F — > F and F — > F — 1 transitions as separate cases. 



1. Dark state for an F — > F transition 

The Clebsch-Gordan coefficients for — F < m < F are given by 



(^1.-^ + 1) = + l ^gf^ )' /2 . (A7) 



Thus, the dark state conditions are 



y/(F + m)(F + 1 - m) a+c m -i = \/(F - m){F + 1 + m) a_c m+ i (A8) 

for — F < m < F, and 

a+CF-i = a-c_(^_!) = 0. (A9) 
For each value of F there is single dark state that meets these conditions. The dark state for a 3 — > 3 transition is 

|d 3 ) =iV(V5a^|3,3) + V3a^a_|3,l> +\/3a + a 2 _|3,-l> +V5a 3 _|3,-3)), (A10) 
and the dark state for a 4 — > 4 transition is 

|d 4 ) = 7V(\/35a+|4,4) + \/20a^a_|4, 2) + \/l8 a+cd |4, 0) + \/20a+ai|4, -2) + V35ai |4, -4», (All) 
where N is a normalization constant that depends on a±. 



Dark states for an F — > F — 1 transition 



The Clcbsch-Gordan coefficients for — F < m < F are given by 



Thus, the dark state conditions are 
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y/(F - m)(F + 1 - to) a + c m -i = -^J{F + m)(F + 1 + to) a_c m+ i. (A14) 

for — F < m < F. For each value of F there are two dark states that meet these conditions. The dark states for a 
4^3 transition are 



\d 4e ) = A^(ct^|4,4) - V28a+a-|4,2) + V70 a 2 + a 2 _\A, 0) - V28 a+a 3 _\4, -2) + ai|4, -4)) (A15) 
|d 4o ) = N(a+|4,3) - v^a+a-l 4 , 1) + V7 a+a^ |4, -1) -ai|4,-3)), (A16) 
where N is a normalization constant that depends on a±. 
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